Particle dynamics and effective temperature of jammed granular 
matter in a slowly sheared 3D Couette cell 
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We report experimental measurements of particle dynamics on slowly sheared granular matter in 
a three-dimensional (3D) Couette cell. A closely-packed ensemble of transparent spherical beads is 
confined by an external pressure and filled with fluid to match both the density and refractive index 
of the beads. This allows us to track tracer particles embedded in the system and obtain three- 
dimensional trajectories, (r(t),9(t), z(t)), as a function of time. We study the PDF of the vertical and 
radial displacements, finding Gaussian and exponential distributions, respectively. For slow shear 
rates, the mean-square fluctuations in all three directions are found to be dependent only on the 
angular displacement of the Couette cell, A0 e , (Az 2 ) ~ A0 e , (Ar 2 ) ~ A0 e a , (AO 2 ) ~ AO/ , where 
a and j3 are constants. With A9 e proportional to the time between measurements, the values of the 
constants, a and /3, are found to be sub-diffusive and super-diffusive, respectively. The linear relation 
between {Az 2 ) and angular displacement implies a diffusive process, from which we can calculate an 
"effective temperature", T e fi, in the vertical direction, through a Fluctuation-Dissipation relation. 
It is of interest to determine whether these systems can be described by analogous equilibrium 
statistical mechanics concepts such as "effective temperature" and "compactivity" . By studying 
the dynamics of tracer particles, we find the effective temperature defined by the Stokes-Einstein 
relation to be independent of the tracer particle characteristic features, such as density and size, 
and dependent only on the packing density of the system. For slow shear rate, both the diffusivity 
and mobility of tracer particles are proportional to the shear rate, giving rise to a constant effective 
temperature, characteristic of the jammed system. We finally discuss the significance of the existence 
of T e g for a statistical mechanics formulation of granular matter. 



I. INTRODUCTION 

Fluctuation-Dissipation (FD) relations are commonly 
used in equilibrium systems, derived from the notion 
that small perturbations and Brownian fluctuations pro- 
duce the same response in a given system [l|. Mobility, 
the constant of proportionality between a particles drift 
speed and a constant external force, is extracted from ve- 
locity statistics of particles in a given system. Diffusivity, 
calculated from fluctuation displacements of particles in 
a system over time, represent the Brownian motion. The 
temperature of a system in thermodynamic equilibrium 
can be extracted from a FD relation, defined as the ra- 
tio of diffusivity and mobility, as is commonly used in 
the Einstein relation. In equilibrium this temperature is 
taken to be the bath temperature. 

As studies in granular matter have grown more im- 
portant within the environmental and industrial fields, 
the need to establish a scientific framework that accu- 
rately predicts granular system responses on the con- 
tinuum level, beyond merely geometrical features, has 
also escalated. Granular matter, when condensed to suf- 
ficiently high volume fractions, undergoes a 'jamming' 
transition to the jammed state. The jammed state is 
defined as the condition when a many-body system is 
blocked in a configuration far from equilibrium, such that 
relaxation cannot occur within a measurable time-scale. 
For granular matter, the jammed state indicates a tran- 
sition between a solid-like behavior, and a liquid-like be- 
havior. At high volume fractions, the physical size of the 
constituent grains inhibits particle motion, thereby ren- 



dering the system out of equilibrium, and the granular 
system behaves more like a solid. Thermal motion does 
not govern the exploration of states in jammed granular 
matter. 

Theories proposed by Edwards and collaborators [3] 
propose a statistical mechanics for granular matter based 
on jamming the constituent grains at a fixed total vol- 
ume such that all microscopic jammed states are equally 
probable and exhibit ergodicity. The exploration of re- 
versible jammed states is achieved via an external pertur- 
bation such as tapping or shear, not Brownian motion as 
in thermal systems. There is an important difference be- 
tween reversible jammed states, and states that are only 
mechanically stable within certain limits of perturbation 
magnitude. For example, pouring grains into a container 
results in a pile at a particular angle of repose. This me- 
chanical equilibrium configuration is jammed regularly 
but not reversibly jammed because in response to an ex- 
ternal perturbation, the constituent particles will irre- 
versibly rearrange, approaching a truly jammed configu- 
ration. Studying an ensemble of truly jammed, reversible 
states is thereby suitable for a plausible application of 
statistical mechanics under the present theory. These 
ensembles, inherently non-equilibrium systems, will not 
be governed by the commonly used parameters of equilib- 
rium statistical mechanics, such as a bath temperature. 

In recent studies theoretical mean-field models of 
glasses [3J have introduced the concept of an "effective 
temperature" as extracted from the FD relations in non- 
equilibrium systems. While not equivalent to the equilib- 
rium bath temperature, the effective temperature reflects 
a change in the relaxation time-scale of the system. These 




FIG. 1: (Color online) Picture of experimental set-up. Trans- 
parent acrylic grains and black tracers in a refractive index 
and density matched solution are confined between the in- 
ner cylinder of radius 5.08cm and the outer cylinder of radius 
6.67cm. 



non-equilibrium systems extend beyond glasses, and into 
granular media, where physical size of the constituent 
grains inhibits motion, allowing for jammed systems far 
from equilibrium. This concept has been furthered by 
computer simulations of granular media and other non- 
equilibrium soft-matter systems [J, d m, 0, @, Q • It re- 
mains a question whether or not granular media can be 
characterized by an effective temperature, thus revealing 
a dynamic counterpart to the static "compactivity" as 
proposed by Edwards @|. 

Athermal systems require the input of energy by an 
external source to explore the effective temperature [l(| ■ 
One proposed method of calculating the effective tem- 
perature of a jammed granula r sy stem is a slow shearing 
procedure [12 0, El El El Efl Ezl Jeading to the de- 
sign of the experiment we present in [18| . Slow shearing, 
at the quasi-static limit, allows for extrapolation towards 
an effective temperature of jammed, static, systems. The 
jammed system of interest is one of identical, spherical 
grains, confined between the two cylindrical walls of a 3D 
Couette cell. The grains are further confined by an exter- 
nal pressure in the vertical direction. The inner cylinder 
of the Couette cell is slowly rotated to induce shearing in 
the system. Tracer particles are inserted in the system, 
and their trajectories recorded via multiple cameras sur- 
rounding the system. The Couette cell is partially filled 
with a refractive index matched fluid to allow for sys- 
tem transparency. The cylindrical walls are roughened 
by gluing grains, identical to those of the bulk, such that 
crystallization is avoided. 

The cameras record tracer particle trajectories 
throughout the bulk, recording data in cylindrical coor- 
dinates, (r(i), 0(i), z(t)). Distributions of tracer particle 
displacements are measured in each direction. As gravity 
is the external force applicable to the mobility calcula- 
tion in the current formalism, only displacements in the 





FIG. 2: (Color online) Top view of experimental set-up. The 
outer cylinder is made of the same material as acrylic grains 
(n ~ 1.49). Once the refractive index is matched, light scat- 
tering from tracers will refract only one times on the outer 
surface of the outer cylinder. A single particle is captured by 
two cameras allowing the determination of the 3-dimensional 
coordinates of the particle, (r,0,z). 



z-direction are applied to the FD relation. Additionally, 
average velocity profiles are calculated for each direction, 
first with constant shear rate, 7 e , and further studied 
to determine shear rate dependence. Displacement mea- 
surements are further limited to the "constant mobility 
and diffusivity" (CMD) region, defined as the narrow 
range of radial coordinates such that the average ver- 
tical velocity is roughly independent on radial distance. 
The PDF of displacement distributions for each direc- 
tion is presented. Further, fluctuations in displacement 
are determined for each cylindrical direction, and stud- 
ied as a function of time. Radial displacement fluctuation 
is found to be sub-diffusive, while angular displacement 
fluctuation is found to be super-diffusive. Vertical dis- 
placement fluctuation is purely diffusive within the time 
scales of the experiment, allowing for the validity of the 
FD relations used herein. All displacement fluctuations 
are reduced to functions of angular displacement, and the 
results are presented. Such relationships permit scaling 
of the PDF curves with varying angular displacements 
due to changing shear rates. 

Utilizing the FD relations presented above, the diffu- 
sivity and mobility in the z-direction are extracted from 
the tracer particle trajectories and the effective temper- 
ature is realized. This effective temperature is found to 
be independent of tracer particle properties, as shown 
in [18j . and further independent of the slow-shear rate. 
Moreover, the effective temperature may then be consid- 
ered a physical variable that characterizes the jammed 




FIG. 3: (Color online) Sketch of experimental set-up. Note 
that the cylinder is surrounded by 4 cameras, in the sketch 
we plot only two cameras. A single particle is captured by 
two cameras allowing the determination of the 3-dimensional 
coordinates of the particle, (r,8,z). 



granular system, with respect to the generalization of 
the equilibrium statistical mechanics of Boltzmann, as 
applied to non-equilibrium systems. 

We further study the limits within which this effec- 
tive temperature may be a valid physical variable, as we 
determine mobility and diffusivity as a function of shear 
rate. While diffusivity appears independent of shear rate, 
even somewhat above the 'slow' regime, mobility shows 
a clear decrease in magnitude as we explore shear rates 
above the slow regime, resulting in an increase in the 
effective temperature as a function of shear rate. 

In this paper, we will further report the experimen- 
tal detail of particle dynamics in [18]. The outline is as 
follows: 



II. EXPERIMENTAL METHOD 

A. Experimental Setup 

The experiment is performed using a three-dimensional 
(3D) Couette cell, as shown in Fig. [U[2]and[3] The grains 
are confined between two cylinders of height 19.0cm. The 
inner cylinder is rotated via a motor, while the outer 
cylinder remains fixed. The walls of the cylinders, in 
contact with the grains, are roughened by means of a 
glued layer of identical granular material, thereby mini- 
mizing wall slip. The walls of the inner and outer cylin- 
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FIG. 4: (Color online) Trajectories of the 3.97mm nylon trac- 
ers in Packing 1 showing the diffusion and response to the 
gravitational force when sheared in the Couette cell. 



ders are roughened by acrylic beads with diameter 3.97 
and 1.59mm, respectively. Testing the experiment with 
a rough inner wall and a smooth outer wall resulted in 
packing crystallization. The grains are compacted by an 
external pressure of a specific value (typically 386 Pa), 
introduced by a moving piston at the top of the granular 
material, acting in the negative z-direction. 

Observation techniques are used to monitor the granu- 
lar packing evolution as it explores the available jammed 
configurations. The Couette cell is sheared at the quasi- 
static limit, with slow frequencies / = 0.2 ~ 4.2 mHz 
defining the external shear-rate j e = 2ir f Ri / (R2 — Ri) = 
uRi/(R 2 -Ri) = 0.004 -0.084 s- 1 , where R t = 5.08cm 
and i?2 = 6.67cm are the radius of the inner and outer 
cylinders, respectively, and u is the angular velocity of 
the inner cylinder (Notice that i?i and R 2 are measured 
after the walls are roughened by a glued layer of beads). 
The experiment is designed to measure the diffusivity 
and mobility of tracer particles 0, LL8|, LL9| , as opposed to 
tracking the motion of all constituent grains. The dis- 
tance between the inner and outer cylinder is less than 
10 g rain diameters to prevent bulk shear band formation 
[12l . Ha , [Fa , [la . [TtJ that may interfere with the experi- 
mental measurements by altering the diffusivity. 

A refractive index matching suspending solution is em- 
ployed in order to create a transparent sample. The sus- 
pending solution is also density matched to the grains in 
order to eliminate pressure gradients derived from gravity 
in the vertical direction, circumventing problems seen in 
previous experiments of compactivity [201 ] and other ef- 
fects such as convection and size seg regation such as the 
Brazil nut effect inside the cell [2l|. The solution used 
in this experiment is approximately 74% weight fraction 
of cyclohexyl bromide and 26% decalin [22(] . These steps 
avoids problems encountered in previous tests of com- 
pactivity. 




FIG. 5: (Color online) A typical trajectory of the 3.97mm 
nylon tracer for 3 hours in 3D plot. The dark gray and light 
gray cylinder indicate the outer surface of the sheared inner 
cylinder and the inner surface of the static outer cylinder re- 
spectively. 
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FIG. 6: (Color online) (a) Average angular velocity, u>e{r), 
(b) Average vertical velocity, v z (r), and (c) Average radial 
velocity, v r (r), versus radial distance r for various tracers 
and different packings. Packing 1 and Packing 2 are run at 
7e = 0.048, 0.024s -1 , respectively. In (a), solid lines are expo- 
nential fitting. In (b), the positive velocity of nylon tracer is 
due to the smaller density than acrylic's. The negative veloc- 
ity of delrin tracer is due to the higher density than acrylic's. 



B. Packing Preparation 
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The granular system is a bidisperse, 1:1 by mass, mix- 
ture of spherical, transparent Poly-methyl methacrylate 
(acrylic) particles, with density p = 1.19 and index of 
refraction n ~ 1.49. The bidisperse mixture is used in an 
effort to inhibit crystallization of the system. The respec- 
tive particle diameters are either 3.17mm and 3.97mm 
(Packing 1) or 3.97mm and 4.76mm (Packing 2). The 
approximate same size ratio of each bidisperse packing 
leads to approximately the same value of volume frac- 
tion for both, being 0.62 before shearing and 0.58 during 
shearing. 

A negative consequence of utilizing a suspending so- 
lution includes possible modification of the friction co- 
efficient between the grains. While this cannot be com- 
pletely avoided, it is important to note that the liquid 
only partially fills the cell (see Fig. [1]), such that the 
pressure of the piston is transmitted to the granular ma- 
terial exclusively, not to the fluid. Additionally, hydro- 
dynamic effects from partial cell filling are avoided by 
the extremely slow rotational speeds applied to our sys- 
tem. The system remains very closely packed, such that 
particles are not free to float in the fluid. Therefore, the 
random motion of the particles is controlled by the 'jam- 
ming' forces exerted by the contacts between neighboring 
grains, not fluid mechanics. 
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FIG. 7: (Color online) Average angular velocity of trac- 
ers, u)$(r), versus radial distance r for various shear rate 
7 C in Packing 2. Black square, red circle, green triangle, 
blue triangle-down, cyan diamond, magenta triangle- left, yel- 
low triangle-right, dark yellow hexagon are corresponding 
to j e = 0.008, 0.016, 0.024, 0.032, 0.041, 0.048, 0.060, 0.084s" 1 , 
respectively. The inset plots the collapsing of average angular 
velocity scaled by shear rate, u$(r)/'j e , versus radial distance 
r for various shear rate j e . The red solid curve is the average 
result of the collapsing. The black solid curve is a exponential 
fitting. 



C. Implementation of Fluctuation-Dissipation 
Theory 

Cylindrical coordinates, (r(t), 9(t), z(t)), of tracer par- 
ticles are obtained by analyzing images acquired by four 
digital cameras surrounding the Couette cell. For sys- 
tems in thermal equilibrium, a Fluctuation-Dissipation 
(FD) relation may be utilized in an effort to calculate 
the bath temperature of the system. This method may 
be extended to non-equilibrium systems, such as jammed 
granular systems presented in this study. The FD rela- 
tion is defined as follows: 



([x(t + At) - x{t)] 2 ) ~ 2D At, 



([x(t + At) - x(t)) ~ MFAt, 



[x(t + At) - x(t)} 2 ) = 2T eff 



(x(t + At) -x(t)) 
F 



(1) 



(2) 



(3) 



The tracer particles must experience a constant force, 
F, in order to calculate the mobility as defined above. 
The most convenient constant, external, force, is gravity 
in the z-direction. If the effective temperature is to be 
regarded as an intensive variable of the non-equilibrium 
system, it requires independence from the tracer parti- 
cles properties, and we present data in favor of this re- 
sult. However, we acknowledge that temperature mea- 
sures from multiple observables would be necessary to 
analyze the underlying thermodynamic meaning of T c g. 



D. Properties of Tracer Particles 

Tracer particles added to the bulk must have proper- 
ties unique from the grains comprising the bulk. How- 
ever, tracer particles too small, or too large, with respect 
to the acrylic grains described previously, would result 
in erroneous measurements. Dynamics of tracer particles 
that were too small would be dominated by "percolation 
effects" [23[ , resulting in larger than expected tracer par- 
ticle displacements. Those too large would require shear 
rates above the quasi-static limit we propose to study, or 
possibly have no dynamics at all due to size limitations. 
With these notions in mind, two different types of tracers, 
nylon (p' = 1.12) and delrin (p' = 1.36), are employed, 
which result in different external forces, F = (p' — p)Vg 7 
where p and p' are the densities of the acrylic particles 
and the tracers, respectively, V is the volume of the tracer 
particle and g the gravitational acceleration. Variations 
in tracer particle diameter and density allow us to study 
dynamical changes due to a change in constant exter- 
nal force, while we remain within a range appropriate to 
achieve results expected to be governed by the effective 
temperature. 
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FIG. 8: (Color online) Average vertical velocity of tracers, 
v z (r), versus radial distance r for various shear rate -y e in 
Packing 2. Black square, red circle, green triangle, blue 
triangle-down, cyan diamond, magenta triangle-left are cor- 
responding to > = 0.008,0.016,0.032,0.048,0.060,0.084s" 1 , 
respectively. The inset plots the collapsing of average vertical 
velocity scaled by shear rate, v z (r)/ A f e , versus radial distance 
r for various shear rate 7 e . The red solid curve is the average 
result of the collapsing. 



E. Particle Tracking Technique 

Four digital cameras symmetrically surround the shear 
cell to track the tracer particles with frame rate ~ 5 
frame/s, as shown in Fig. [3] The outer cylinder is made 
of the same material as the grains (acrylic, n ~ 1.49). 
The refractive index is matched by the fluid such that 
the system can be regarded as an optical whole, i.e., the 
light scattered from tracers refracts only once at the outer 
surface of the outer cylinder, as shown in Fig. [5] The 
determination of the 3D tracer position is achieved by a 
simple calculation considering both system geometry and 
2D projections captured by two adjacent cameras. 

Camera calibration and determination of relative po- 
sition is important as a minimal asymmetry will result 
in a large calculation error of the tracer particles coordi- 
nates, (r(t), 8(t), z(t)). As opposed to directly measuring 
relative positions of the cameras by physical devices, we 
utilize computer programming. In order to simplify the 
calculation, we assume the camera to be a pinhole, mean- 
ing all light coming into the camera coincides at a single 
focal point. 

Before each experiment, we record the images of a piece 
of grid paper attached to the surface of the outer cylin- 
der, acting as the 2D projections of the outer cylinder for 
each camera. Next, we adjust the positions of four cam- 
eras until each camera can give the approximately same 
2D projections of the outer cylinder. Then we use the 
computer program to generate a virtual cylinder, along 
with four virtual cameras, according to the geometry of 
the shear cell, In other words, we build a virtual space 
of the entire experimental setup and the respective geo- 



metrical relations between its elements. 

From the previous calibration procedure, we have the 
relative positions of the four cameras to the shear cylinder 
with sufficient accuracy. In order to further calibrate and 
know the exact position of cameras, we adjust the relative 
position of cameras in our virtual space until the virtual 
2D projection of the cylinder to the cameras coincides 
exactly with the actual projection, being the grid paper 
attached to the outer cylinder. When this procedure is 
accomplished, the virtual space exactly coincides with 
the real experimental setup space. Therefore the virtual 
relative position of cameras are also the real positions. 

Furthermore, in our virtual space, any point with 3D 
coordinates, (r, 8, z), we can calculate its 2D coordinates 
in four virtual 2D projections, {x\, y\) ~ (#4, 2/4), by con- 
sidering the geometry relation to cameras. Oppositely, 
for any tracer particle, if we know its 2D coordinates in 
four 2D projections, (xi,y{) ~ (£4,3/4), we can exactly 
locate its 3D positions, (r, 0, z), since the virtual space 
is equal to the actual one. The resulting vertical trajec- 
tories of the tracers z(t) are depicted in Fig. [1] showing 
that the nylon tracers not only diffuse, but also move 
with a constant average velocity to the top of the cell. 
Fig. [5] shows a typical trajectory of tracer particle in 3D 
plotting. 
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FIG. 9: (Color online) Average radial velocity of tracers, 
v r (r), versus radial distance r for various shear rate 7 e in 
Packing 2. Black square, red circle, green triangle, blue 
triangle-down, cyan diamond, magenta triangle- left, yellow 
triangle-right, dark yellow hexagon are corresponding to j e = 
0.008, 0.016, 0.024, 0.032, 0.041, 0.048, 0.060, 0.084s" 1 , respec- 
tively. The inset plots the collapsing of average radial veloc- 
ity scaled by shear rate, w r (r)/7 e , versus radial distance r for 
various shear rate y e . The red solid curve is the average result 
of the collapsing. 



III. RESULTS 
A. Average Velocity Profiles 

We first study the velocity profiles for a fixed shear 
rate, j e , followed by a study on the shear rate depen- 
dence in the next section. The average velocity profiles 
in the angular direction, ujg(r), in the vertical direction, 
v z (r), and in the radial direction, v r (r), are obtained 
by averaging the velocities of all tracer particles over all 
times at each radius r, as shown in Fig. [5] 

As observed in previous work [14|, we find that uig{r) 
can be expressed in the exponential form demonstrated 
in Fig. [5^,: 



u e (r) = Ai 



7e(i?2 



Rl 



Ri) f . r 
exp(-A 2 



Ri 



i?2 - R\ 



Aiwexp(— A2 



Ri 



(4) 



i?2 — Ri 



where Ai and A2 are constants independent of shear rate, 
tracer size and tracer type, depending only on the type 
of packings and geometry of the shear cell. We find 
that Ai = 0.77 and 0.73, and A 2 = 2.15 and 1.43 for 
Packing 1 and Packing 2, respectively. When r = R\, 
we(Ri) = Aiw being the angular velocity of the first layer 
of grains closest to the sheared inner cylinder with angu- 
lar velocity 10. Therefore Ai (0 < Ai < 1) can be taken as 
the efficiency of shearing, describing the amount of slip 
between the inner rotating cylinder and the first layer of 



grains it contacts. Packing 1 has a higher value of Ai 
than Packing 2, as the smaller grains follow the rotat- 
ing inner cylinder more easily. On the other hand, when 
r = i?2, we find u>g{R2) — Ai^exp(— A2), the velocity of 
the last layer of grains closest to the static outer cylin- 
der. This velocity is non-zero, so that the shear band is 
located right at the outer cylinder, avoiding the forma- 
tion of shear bands in the bulk. 

In order to avoid the tracer particles sticking to the 
outer cylinder surface and forcing its velocity to zero, we 
glue smaller size particles to roughen the outer cylinder. 
This roughens the surface of the outer cylinder and avoids 
crystallization. Further, this allows slipping of the bulk 
particles at the outer cylinder, forcing the shear band to 
be located exactly at the outer cylinder, not in the bulk. 
The glued particles are 1.59mm, smaller than the sheared 
granular material. The mean angular and vertical veloc- 
ity, uig and v z , of tracer particles do not decay to zero 
even if the tracers come close to outer cylinder surface. 
(See Fig. [7]and Fig. [S]at r = R 2 ). 

The exponential decay of ujg(r) results in a local shear 



rate, 



diAjg(r) 



,, , dependent on radial distance. In the re- 
gion near the outer cylinder, U)$(r) decays slowly with 
increasing r which leads to weak dependence of "i 
on the radial distance r. If we do Taylor expansion 
at r = i?2) the average angular velocity of the tracers, 
ujg{r), can be approximated to a linear function of r, 



u}g(r) fa Aiwexp(-A 2 ) - r 



dr 



with constant local 



shear rate ^^ = 0.021 s^cm" 1 



The diffusivity and 
mobility of the tracer particles strongly depend on the 
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FIG. 10: (Color online) (a) PDF of the vertical displacements, 
P(Az), of the 3.17mm delrin tracers in Packing 1 for a given 
time interval At = 50s, and with 7 e = 0.048s -1 . Tracer tra- 
jectories are split into sub-trajectories confined in two regions, 
(i): 5.08cm < r < 5.80cm, which is close to inner rotating 
cylinder, and (ii): 5.80cm < r < 6.67cm, which is far away 
from inner rotation cylinder. We compared the calculated 
P(Az) by using the sub-trajectories from the regions of (i) 
and (ii) respectively, which are plotted as black triangle and 
black circle. See more details in the main text, (b) PDF of the 
vertical displacements, P(Az), of the 3.97mm nylon tracers 
in Packing 1 with j e — 0.048s -1 , shifted by the average dis- 
placement (Az) and scaled by the root-mean-square deviation 
(Az(t) 2 ) 1 ' 2 . The red solid curve is a Gaussian distribution, 
P{x) = 0Ae~ x2/2 . 



local rearrangement of the grains. A constant shear rate 
results in homogenous local rearrangement of the pack- 
ings ensuring that the diffusivity and mobility of tracers, 
dependent on local shear rate, remain approximately in- 
dependent of r. As shown in Fig. [5b, we find a plateau 
in the vertical velocity profile which can be further seen 
in Fig. [5J Similar behavior is observed in the vertical 
diffusivity profile, D z (r) as shown in Fig. 1 161 which we 
will discuss in detail in Section IIIIDI We denote this 
the "constant mobility and diffusivity region", i.e., CMD 
region, 5.80cm < r < 6.67cm. Contrary to prior work 
pj], El E! El ES El E3 on sheared granular matter in 
the Couette cell, our experiment focuses only a narrow 
gap, 15.9mm, of Couette cell. The CMD region allows 
us to well define the diffusivity and mobility of the tracer 
particles, such that we can calculate the average vertical 
velocity, v Zl and the average vertical diffusivity, D z , by 
averaging the velocities of all tracers over all times in the 
CMD region, significantly improving the statistics. In 
this study, the statistical average and the measurements 
of tracer fluctuations will be confined only to the CMD 
region. 

We find v r (r) to be flat for different types of tracer 
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FIG. 11: (Color online) PDF of the radial displacements, 
P(Ar), of the 3.97mm nylon tracers in Packing 1 with 7 e = 
0.048s -1 for different time intervals. A symmetric distribu- 
tion around zero displacement indicates that there is no net 
flow in the radial direction. The solid lines are exponential 
fitting, P{Ar) ~ exp(— l^i), where r = 0.17,0.25,0.32 for 
At = 100, 300, 500s, respectively. The inset shows the rms 
fluctuations, which gives the value of a = 0.67. 



particles and for different packings except when the trac- 
ers are close to the inner and outer cylinder, i.e., r = R\ 
and r = i?2, as shown in Fig. [6fc. V r (Ri) is negative and 
v r {R2) is positive, indicating the inner and outer cylin- 
der walls can slightly attract the tracers, ft should be 
noted that the statistics presented in this study does not 
incorporate data from the regions close to the inner and 
outer cylinder to avoid these boundary effects. 



Shear Rate Dependent Average Velocity 
Profiles 



Next, we study the dependence of the particle velocity 
on the external shear rate. According to Eq. (0|, the 
velocity profile in the angular direction, u)g(r), is pro- 
portional to the external shear rate j e . We can collapse 
to$(r) by scaling the shear rate. The results are shown 
in the inset of Fig. [7] for Packing 2. The collapsing of 
ujg(r)/j e shows a periodic shape superimposed to expo- 
nential decay with a very small amplitude, also found in 
the velocity profile of v r (r) (see Fig. ^. The periodic 
length is roughly equal to the grain particle size and re- 
flects the different layers of grains in the radial direction. 
This periodicity is weaker in Packing 1 than Packing 2, 
since the particle size of Packing 1 is smaller than that 
of Packing 2. 

The collapsing method can be further applied to v z (r), 
as seen in Fig. [8] After scaling by the shear rate, v z (r)/j e 
also shows a flat plateau indicating the CMD region. 
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FIG. 12: (Color online) PDF of the angular displacements, 
P(A6), of the 3.97mm nylon tracers in Packing 1 with y e — 
0.048s _ for different time intervals. Due to Taylor dispersion 
effects the distribution shows an asymmetric shape. The rms 
fluctuations shown in the inset reveal a faster than diffusion 
process. 



C. Probability Distribution of Displacements 
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FIG. 13: (Color online) PDF of the vertical displacements, 
P(Az), of the 4.76mm nylon tracers in Packing 2 for various 
effective angular displacement A9 e and effective shear rate 7 e . 
The PDFs are scaled by AOe 1 ' 2 and shifted by the mean dis- 
placement (Az). The red dashed line is the Gaussian fitting, 
P(x) ~ exp[— ( ^ 47 ) 2 ]. The collapsing of PDFs indicates that 
the RMS fluctuations of the vertical displacements follow the 
relation, (Az 2 ) ~ A0 e - 



Fig. [10b shows the results of the probability distribu- 
tion of the displacements Az in the vertical direction for 
a given time interval At. The data corresponds to the 
3.17mm delrin tracers in Packing 1. Usually, 20 tracers 
are used for calculations. Tracer trajectories are split 
into sub-trajectories confined in two regions, (i): 5.08cm 
< r < 5.80cm, close to the inner rotating cylinder, and 
(ii): 5.80cm < r < 6.67cm, i.e., CMD region, close to 
the outer cylinder. We compare the calculated P(Az) 
by using the sub-trajectories from the regions of (i) and 
(ii) respectively, which are plotted as black triangle and 
black circle in the Fig. [10k . The data in the inner region 
(i) clearly display an asymmetric tail for Az < 0. This 
extra spreading is similar to the phenomena of the Taylor 
dispersion [24J . 

Taylor dispersion appears when diffusion couples with 
the gradient of flow giving rise to a larger dispersion along 
the flowing direction (see for instance [T7j for a study of 
Taylor dispersion in granular materials). In the present 
experiment, the shear rate of granular flow in the angular 
direction exhibits exponential decay, as shown in Fig. 
O). The larger shear rate in the inner region (i) results 
in larger packing rearrangement, which gives rise to a 
larger dispersion in the vertical direction. In this case 
it is not possible to extract the bare diffusion constant. 
On the contrary, for the region (ii), i.e., CMD region, as 
we mentioned, the gradient of the flow, i.e., the shear 
rate is approximately constant, giving rise to a Gaussian 
diffusion, as shown in Fig. 110a . By measuring the width 
and the mean value of this Gaussian distribution of 
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exp[ 



exp[ 



(Az-(Az)) 2 
2<Az 2 } J 
(Az-M z FAt) 2 
4D?At ' 



(5) 



we can define the diffusivity and mobility, D z and M Z1 
which lead to the effective temperature of the granular 
packing discussed in the following section. In the Fig. 



[TOb . we define a new scaled variable 



__ Az-(Az) 



and 



(A Z 2>!/2 

plot P(x) for different At, all the curves are found to 
collapse into a single curve 



P(x) 



3 -x 2 /2 



(6) 



In this experiment, we will focus our measurements in 
the region away from the inner boundary (region (ii), i.e., 
CMD region), where the mobility is a constant (as shown 
a plateau in the inset of Fig. [5]) and Taylor dispersion 
effects are absent. 

We find exponential fluctuations for the probability 
distributions of the tracer particles in the radial direc- 
tion as shown in Fig. [Ill 



P(Ar) 



(7) 



where r is a function of At. The symmetric shape for 
P(Ar) indicates the absence of a shear induced segrega- 
tion, as observed with multiple sizes of grains, as there 




FIG. 14: (Color online) PDF of the radial displacements, 
P(Ar), of the 4.76mm nylon tracers in Packing 2 for various 
effective angular displacement A9 e and effective shear rate j e . 
The PDFs are scaled by A6 e a/2 , where a = 0.67. The red 
dashed line is exponential fitting, P(x) ~ exp(— _ q1 9 ). The 
collapsing of PDFs indicates that the RMS fluctuations of the 
radial displacements follow the relation, (Ar 2 ) ~ A8 e a . 



is no net flow of the tracer particles towards either cylin- 
drical wall within the time-scales of the experiment. We 
also observe no average motion of the tracer particles 
towards the center of the Couette cell except within a 
small range of radial distance, around 0.12cm, close to 
both walls where particles experience a slight attraction 
to boundaries. These features are shown in Fig. [6b and 
Fig. m 

The analysis of the radial displacement fluctuation re- 
veals a power law, sub-diffusive, process: 



(Ar 2 ) ~ At a 



(8) 



as shown in Figure 8, where a — 0.67 for both the delrin 
and nylon tracers. 

The data taken for the angular displacement is in the 
direction of the flow, and affected by Taylor dispersion 
as shown in the non-Gaussian tail of the displacement 
distribution A6(t) in Fig. Q21 This leads to a power law, 
super-diffusive process, illustrated by 



(Ae 2 ) ~ At' 3 



(9) 



as seen in the analysis of the fluctuations of A8 shown 
in Fig. [121 where j3 = 1.67 and 1.30 for Packing 1 and 
Packing 2 respectively. 

We further study how shear rate affects the displace- 
ment probability distribution. We find that for small 
shear rate, the probability distributions of displacements 
in the three cylindrical coordinates are independent of 
the shear rate, depending only on the sheared displace- 
ment, i.e., the external rotating displacement, defined as 
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X = (A0-<A8>)/A8 



FIG. 15: (Color online) PDF of the angular displacements, 
P(A6), of the 4.76mm nylon tracers in Packing 2 for various 
effective angular displacement A9 e and effective shear rate 
7 e . The PDFs are shifted by the mean displacement (Ad) 
and scaled by AO/' 2 , where f3 — 1.30. The red dashed line 
are Gaussian and exponential fittings for x < and x > 0, 
respectively. The collapsing of PDFs indicates that the rms 
fluctuations of the angular displacements follow the relation, 
(A<9 2 ) ~ A9/. 



A9 e = j e At 

= u:AtR 1 /{R 2 - Ri) 
= A9 l R 1 /(R 2 - R t ) 



(10) 



where A6i is the rotating displacement of the inner cylin- 
der. This result is expected. Since we shear the Couette 
cell very slowly, the diffusion of the tracers depends only 
the number of granular packing configurations sampled 
by the Couette cell, which depends only on the sheared 
displacement. 

As emphasized in the previous text, the statistical av- 
erage and the measurements of the tracer fluctuations is 
confined to the CMD region, such as the D z shown in Fig. 
[PTk . We calculate the D z by measuring the width and the 
mean value of the Gaussian distribution of P(Az), and 
obtain the P(Az) by averaging the displacement fluctua- 
tions of all tracers over all time in the CMD region. Next, 
we apply a different method to reveal how D z (r) depends 
on the radial distance r, as shown in the Fig. [TCJ We 
first obtain P(Az, r) for a certain radial distance r, then 
we calculate D z (r) by measuring the width and the mean 
value of the Gaussian distribution of P(Az,r). In Fig. 
[Pol we see that the tracer particles have higher diffusivity 
close to the inner cylinder than the outer. Since D z ~ j e , 
we can collapse all the D z [r) for various shear rates, as 
shown in the inset of Fig. [121 The collapse of D z (r)/% 
shows a plateau close to the outer cylinder, consistent 
with our previous discussion of the CMD region. 
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FIG. 16: (Color online) Diffusivity D z versus radial dis- 
tance r for various shear rate 7 e in Packing 2. Black 
square, red circle, green triangle, blue triangle-down, cyan 
diamond, magenta triangle-left are corresponding to 7 e = 
0.008, 0.016, 0.032, 0.048, 0.060, 0.084s -1 , respectively. The 
inset plots the collapsing of diffusivity scaled by shear rate, 
Dz/^e, versus radial distance r for various shear rate j e . The 
red solid curve is the average result of the collapsing. 
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Eq. (fTTj) implies that one can collapse the probability dis- 
tribution of the displacements, P(Az), P(Ar) and P(A6) 
for different shear rates and time intervals by scaling Az, 
Ar and A6 respectively to Az/A6 e 1/2 , Ar/A8 e a/2 and 
A6/A6 e p/2 . The results are presented in Fig. [TO] Poland 

m 



D. Effective Temperature 

We present results for the diffusivity in the z direction, 
the only direction where the effective temperature can be 
calculated due to the vertically acting external force. The 
Gaussian distribution in P(Az) allows us to apply the FD 
relation to the particle displacements, as the diffusivity is 
proportional to the variance of a Gaussian distribution in 
displacements. Exponential fluctuations do not possess 
this same property, but it is important to note that the 
radial direction has no constant applied external force. It 
remains a possibility that a well-defined effective temper- 
ature for displacements in the radial direction could exist. 
To test whether the effective temperature is isotropic, as 
done in [25j| . may be of great interest in future studies. 

A common method of performing a time average to 
measure transport coefficients is employed (see Chapter 
5.3 in [26() by dividing the trajectory of a single tracer 
particle into a series of trajectories, having evenly spaced 



start times, separated by time interval At. The diffu- 
sion constant is obtained by averaging over the aggre- 
gate of tracers and over the initial time intervals, allow- 
ing for the use of merely 20 tracer particles in this par- 
ticular system. Correlations between measurements are 
ensured to have decayed almost to zero, rendering time- 
translational invariance valid in this system, without any 
measurable "aging", since under shearing, system reaches 
the "stationary state" J27J. Furthermore, doubling the 
number of tracer particles leaves D z unchanged, indicat- 
ing independence of the diffusion constant from the num- 
ber of tracers that explore the jammed configurations of 
this non-equilibrium system. 

Analysis of the vertical particle displacements in the 
CMD region reveals a Gaussian distribution, broadening 
over time, as seen in Fig. [TO] and Fig. [TO] For suffi- 
ciently long times period, the mean square fluctuations 
grow linearly (see Fig. [PTk): 



([z(t + At) - z{t)f) ~ 2D z At, 



(12) 



where D z is the self-diffusion constant in the vertical di- 
rection. For the both nylon and delrin 3.97mm tracers 
in Packing 1 we obtain D z 3.97 mm ~ (1.15 ± 0.1) x 10~ 8 
m 2 /s. 

Figure [TTb shows mean value tracer particle positions, 
extracted from the peak of the Gaussian distribution, as 
a function of time. The nrobility in the vertical direction, 
M z , is defined as 



(z(t + At) - z(t)) ~ M z FAt. 



(13) 



The applied force on the tracers, F = (p — p')Vg, is the 
gravitational force due to density mismatch where p and 
p' are the densities of the acrylic particles and the tracers, 
respectively, V is the volume of the tracer particle and g 
the gravitational acceleration. The value of the mobility 
for the both nylon and delrin 3.97mm tracers in Packing 

1 is M z 3.97mm » (9.4 ± 0.9) X 10" 2 s/kg. 

Fig. 117b further reveals a downward curvature of the 
mean-square fluctuations, for sufficiently long times pe- 
riod. Additionally, an apparent cut-off time for the tracer 
particles fluctuation measurements is shown. These ef- 
fects are due to the finite size effect imposed upon the 
tracers by the finite trajectories and should be inversely 
proportional to the tracer particles velocities. Tracer par- 
ticles with larger mobility will have larger mean velocities 
and take a shorter time to complete its trajectory in the 
cell. The cut-off discussed in reference to Fig. [TTk is 
prominently displayed in the 3.17mm delrin tracers of 
Packing 1, having the largest mobility, hence increased 
mean velocities, as shown in Fig. 117b . The larger mobil- 
ity results in the shortest cut-off time for the diffusivity. 
Conversely, 3.97mm delrin tracers of Packing 1 have a 
smaller mobility, hence a longer cut-off time for the dif- 
fusivity. It is important to note that for all tracer par- 
ticles studied here, the cut-off is observed for distances 
larger than a few particles diameters, ensuring that the 
study examines the structural motion of the grains and 
not internal nrotion inside of "cages" . 
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Packing 1, 1:1 mass mixture of 3.17mm & 3.97mm acrylic beads, 7 e = 0.048s -1 


tracer 


d (mm) 


P (g-cm" 3 ) 


D z (10- 8 m 2 *- 1 ) 


M z (io- 2 s-kg- 1 ) 


Tcff (io- 7 J) 


acrylic 


3.17 


1.19 


2.5 ±0.3 






delrin 


3.17 


1.36 


2.4 ±0.3 


24 ±3 


1.0 ±0.2 


delrin 


3.97 


1.36 


1.2 ±0.1 


9.3 ±0.9 


1.3 ±0.2 


nylon 


3.97 


1.12 


1.1 ±0.1 


9.5 ±0.9 


1.2 ±0.2 


ceramic 


3.97 


3.28 




2.2 ±0.2 




brass 


3.97 


8.4 




1.7 ±0.1 




Packing 2, 1:1 mass mixture of 3.97mm & 4.76mm acrylic beads, ■je = 0.024s" 1 


tracer 


u (mm) 


P (g-cm 3 ) 


D z (io- s m 2 *- 1 ) 


M z (io- 2 s-kg- 1 ) 


T cff (io- 7 j) 


nylon 


3.97 


1.12 


1.8 ±0.1 


19.0 ±0.9 


0.95 ±0.07 


nylon 


4.76 


1.12 


1.6 ±0.1 


15.7 ±0.4 


1.0 ±0.1 



TABLE I: Diffusivity and mobility for the different types of tracer and packings. 
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<] 3.17mm acrylic (packingl) V 3.97mm nylon (packing) 

D 3.17mm delrin (packingl) 1.7(>mm nylon ( packing 1 

O 3.97mm nylon (packingl) 

±- 3.97mm delrin (packingl) 




□ 3.17mm delrin (packingl ) 
O 3.97mm nylon (packingl) 

3.97mm delrin (packingl) 
V 3.97mm nylon (packing2) 

4.76mm nylon (packing2) 
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\7 3.97mm nylon (packing2) 






4.76mm nylon (packing2) 
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responses, with At, as the parameter, as extracted from 
Fig. [TTr and Fig. [TTb . A linear relationship exists be- 
tween diffusivity and mobility, with a slope of T c g. We 
obtain for the both nylon and delrin 3.97mm tracers in 
Packing 1, T eff w (1.25 ± 0.2) x 1CT 7 J. 

If the effective temperature is to be regarded as an 
intensive thermodynamic quantity, changing the tracer 
particle size should give rise to a different diffusion and 
mobility yet result in the same measurement of effective 
temperature. The above calculation is repeated for delrin 
tracers of 3.17mm in Packing 1. We find that while the 
mobility and diffusivity change dramatically with respect 
to tracers of 3.97mm, (D z 3 .i 7mm = (2.4 ± 0.3) x IO -8 
m 2 /s and M z 3.i7 mm = (2.4 ±0.3) x 10 _1 s/kg) as shown 
in Table HI due to the change in tracer size, their ratio 
remains unchanged. In all cases D z and M z are inversely 
proportional to the size of the tracers, but the effective 
temperature remains approximately the same, as seen in 
Fig. [TTb . with an average value over all tracers of 



T cff « (1.1 ±0.2) x 10" 'J 



(15) 



FIG. 17: (Color online) (a) Autocorrelation function of trac- 
ers, (b) Response function of tracers, (c) Log-log plot of 
effective temperatures for various tracers and different pack- 
ings as obtained from a parametric plot of their autocorrela- 
tion function versus response function, (d) Same as (c) but 
in a linear-linear plot. The slopes for different tracer diffu- 
sivity vs. mobility curves return the same average value of 
T eS « (1.1 ± 0.2) x 10~ 7 J as given by Eq. ©. 



According to a Fluctuation-Dissipation relation, we 
calculate T e g: 

F{[z{t + M)-z{t)f) 
cff 2(z(t + At) - z(t)) ' [ ' 

Fig. [TTb shows a parametric plot of fluctuations and 



Though this effective temperature is high with respect 
to the bath temperature, we note that a plausible scale 
for the system energy [1Q|, is (p — p')gd, the gravita- 
tional potential energy to move a nylon tracer parti- 
cle one particle diameter, d. A corresponding temper- 
ature would arise from the conversion of this energy 
into a temperature via the Boltzmann constant, ks, is 
T cff = 2.7 x IO 13 k B T at room temperature (T = 300K). 
This specific value serves as a coarse-grained estimate, 
since the tracer size and density clearly shift its value, 
and we focus on the order of magnitude. This large value 
is expected [10], and agrees with computer simulation es- 
timates for an athermal granular system [4|. Therefore, 
our calculated value for T e ff in a sheared granular system 
appears reasonable within the boundaries of the present 
theory. 
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FIG. 18: (Color online) Mobility, M z , versus the external 
force, F, for 3.97mm tracers in Packing 1. The black squares 
are experimental data which is coming from different types 
of tracer, they are (from left to right) nylon, delrin, ce- 
ramic and brass. The red dashed guide line is a function 
of 8.0exp[-(x/0.007) 4 ] + 1.7. 



E. Linear Response Regime 

In an effort to further test the concept of the effective 
temperature as an intensive quantity, a linear response 
regime in the system is of great interest. Such a regime 
would imply that mobility and diffusivity are indepen- 
dent of the external gravitational force as F — ► 0. The 
external force is varied by changing the density of the 
tracers of the same size. This is realized experimentally 
in Packing 1 by the introduction of delrin (p' = 1.36) 
tracers of 3.97mm diameter, the density of which is higher 
than that of nylon (// = 1.12). 

Analysis of the trajectories reveals that the mobility 
is approximately the same for both the delrin and nylon 
tracers with the same diameter and is thereby indepen- 
dent of the external force, as shown in Fig. [TTb and Fig. 
IT51 Further, the external force should have no effect on 
the diffusivity. By calculating the diffusivity of the non- 
tracer particles via dying acrylic tracers and analyzing 
their trajectories, as shown in Fig. [TTb . the diffusion of 
the acrylic tracers of size 3.17mm (for which no external 
force is applied) is the same as the diffusion of the delrin 
tracers of the same size (for which the gravitational force 
is applied). A further example of this property would be 
using two different tracers with different sizes, but hav- 
ing the same external force applied. One would calculate 
two different values of the diffusivity, due to the variation 
in tracer size, without having any variation in external 
force. 

Nonlinear effects appear for tracers heavier than del- 
rin, implying that mobility depends on the external force 
for large enough forces. We find that for a 3.97mm ce- 
ramic tracer (// = 3.28) in Packing 1 the mobility is 



Mz ceramic = (2.2 ±0.2) x 10 2 s/kg and for a brass tracer 
(// = 8.4), M z brass = (1-7 ± 0.1) x 10~ 2 s/kg as shown 
in Table |TJ smaller than the mobility of the nylon and 
delrin tracers of the same size. This behavior is expected 
since if a linear regime exists in the system, it will be 
valid only within certain limits, i.e., M z remains a con- 
stant for small value of external force, F, as shown in Fig. 
1181 It is here that our experiments approach the bound- 
aries presented above for estimated of energy scales for 
the sheared granular system. Our effective temperature 
measurements are therefore limited to those tracer parti- 
cles for which we experience a linear regime with respect 
to both mobility and diffusivity. 

Lastly, the experiment is again repeated for a different 
packing of spherical particles, noted earlier as Packing 
2. Having nearly the same volume fraction of particles 
being both packings of spherical particles, one would ex- 
pect T e g to remain unchanged, as it is a measure of how 
dense the particulate packing is (i.e. a large T c q implies 
a loose configuration, e.g. random loose packing, while a 
reduced T e g implies a more compact structure, e.g. ran- 
dom close packing) . It is found that although differences 
exist between the two packings with respect to mobility 
and diffusivity, as shown in Fig. [T7b.b. the effective tem- 
perature remains approximately the same, as shown in 
Fig. [TTb . It should be noted that both packings are com- 
posed of spherical particles and the statement regarding 
effective temperature as a measure of particulate packing 
density would not be true if the packings are composed 
of particles of, for instance, different shapes, even if they 
have the same volume fraction. 



F. Shear Rate Dependence 

We further explore the effective temperature as an in- 
tensive quantity by analyzing diffusivity and mobility as 
a function of the shear rate. We show in Fig. [TIJ] that 
the effective temperature seems to become approximately 
constant, as long as the particulate motion is slow enough 
such that the system is very close to jamming. We find 
that 



D, 



M z 



7 e < 0.06s" 



7e 



7. 



< 0.04s" 



(16) 
(17) 



while T c ff = D z /M z remains approximately constant for 
sufficiently small j e . 

It is within this quasi-static range where the effective 
temperature could be identified with exploration of the 
jammed configurations. As it remains an important as- 
sumption of this study that the system is being contin- 
uously jammed, shear rates high enough to impact the 
effective temperature measurement imply systems that 
are not continuously exploring jammed configurations. 
As we study the nature of the jammed granular pack- 
ings, it is logical to presume that quasi-static shearing 
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FIG. 19: The dependence of (a) diffusivity D z , (b) mobility 
M z and (c) effective temperature T e g on the shear rate -y e 
for the 4.76mm nylon tracers in Packing 2. The solid lines 
in (a) and (b) are linear fitting [note that the line in (b) is 
a fitting only for the first 6 data points at the small value of 
shear rate j e ]. We find that D z ~ 7 e and M z ~ j e , while 
T c ff = D z /M z is approximately constant for sufficiently small 
j e . This quasi-static regime coincides with the appearance 
of a rate-independent stress in experiments [29J ]. that T e g- is 
interpreted as the temperature of the jammed states. The 
height of flat solid line in (c) is calculated from the slope 
of lines in (a) and (b), which indicates a constant effective 



temperature T e g 
shear rate j e . 



(1.2 ± 0.2) x lCT'J at the small value of 



will provide systems of interest. The limit of T c g as 
7 e — > may result in an effective temperature for the 
static jammed configuration. The quasi-static shear rate 
regime observed could be analogous to the shear-rate in- 
dependent regime observed in the behavior of shear stress 
in slowly sheared granular materials [28|, [29( . This solid 
friction- like behavior has been previously studied [28|, [29| 
and occurs when frictional forces and enduring contacts 
dominate the dynamics. This regime has been also ob- 
served in recent computer simulations of the effective 
temperature of sheared granular materials @, [3CJ . Our 
calculations of T c g for systems close to jamming exclude 
the systems outside of the quasi-static range, in accor- 
dance with prior studies. 



IV. OUTLOOK: SIGNIFICANCE OF T eff FOR A 
STATISTICAL MECHANICS OF GRAINS 

In contrast to measurements of slow mode tempera- 
tures, exemplified by T e ff, we also measure the tempera- 
ture of the fast modes as given by the root mean square 
(RMS) fluctuations of the velocity of the particles. It 
should be noted that these velocities are not instanta- 
neous, as the time necessary to obtain an instantaneous 



velocity is much smaller than the time between measure- 
ments. Nevertheless, we can obtain an estimate of the 
kinetic granular temperature, Tk, from T k = §£*;, where 
Ek = \rriv 2 with v 2 the average kinetic energy of the 
grains. We obtain T k = 9.17 X 10 10 fc B T, or 3.77 x 1(T 10 J 
and T k = 1.54 x 10 n fc B T, or 6.34 x 1CT 10 J, for 3.17 
mm and 3.97 mm delrin tracers in packing 1, respec- 
tively. Here, T = 298.15iC, the room temperature, and 
k B = 1.3806504 x 10~ 23 JK~ X . This kinetic granular 
temperature is smaller than T e g and differs for each type 
of tracer indicating that it is not governed by the same 
statistics. Similar results have been obtained in exper- 
iments of vibrated granular gases [311 ] . The significance 
of this result is that fast modes of relaxation are gov- 
erned by a different temperature. This result is analo- 
gous to what is found in models of glasses and computer 
simulations of molecular glasses (see for instance Refs. 
[1 H3, H [H ii|). In the glassy phase of these models, 
the bath temperature is found to control the fast modes 
of relaxation and a different, larger, effective temperature 
is found to control the slow modes of relaxation. Simi- 
larly, we find a granular bath temperature for the fast 
modes and a larger effective temperature for slow modes 
of relaxation. 

It is possible to identify T e ff as the property of the 
system governing the exploration of jammed configura- 
tions. As this particular non-equilibrium system remains 
athermal, the 'bath' temperature in which the grains ex- 
ist is immaterial, as shown above. Particle diffusion is 
of the order of several particle diameters over the time 
scale of the experiment (see Fig. 2] and [TTa ) implying 
that exploration of the available jammed configurations 
occurs via rearrangements of the particles outside their 
"cages" , suggesting that the trajectory of the system can 
be mapped onto successive jammed configurations ex- 
plored by the system. 

Incorporating certain experimental conditions of re- 
versibility, and ergodicity, a statistical mechanics for- 
mulation may well describe a jammed granular system 
0, (35|. Under the primary assumption that different 
jammed configurations are taken to have equal statistical 
weight, observables can be calculated by "flat" averages 
over the jammed configurational space [2. |4J. 171 13a. I37L l38l. 
[39| . This assumption, advocated by Edwards and collab- 
orators, has been thoroughly debated in the literature 
(see for instance [271 . |4Q | ) . Existing work suggests the ef- 
fective temperature obtained by applying a fluctuation- 
dissipation theory to non-equilibrium systems is analo- 
gous to performing a "flat" average over the jammed 
configurational space, at least for frictionless systems J4J. 
Additionally, the effective temperature can be identified 
with the compactivity introduced in [2j, resulting from 
entropic calculations of the granular packing [J, @, l38| . 
Experimentally testing these ideas is difficult as the en- 
tropy of the jammed configurations is not easily mea- 
sured, and it is not possible to obtain the compactivity 
from entropic considerations in the present study. 

The exploration of reversible jammed states in gran- 
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ular matter bears similarity to that of inherent struc- 
tures in glasses. Inherent structures form a network of 
attractive basins within an energy landscape, and the 
system explores these basins as governed by their stabil- 
ity over the slow-relaxation time of the glass. It should 
be noted, however, that there exists a crucial difference 
between glasses and grains. In liquids energy remains 
conserved, while energy is dissipated in granular systems 
through frictional contact and path dependent forces be- 
tween grains. Thus, a driven granular system will quickly 
come to a mechanically stable, or jammed, state after the 
removal of the driving forces. By its nature, energy is not 
conserved in a granular system. As energy conservation 
is the crucial property used to define an energy ensem- 
ble in statistical mechanics, the use of energy to char- 
acterize granular systems is questionable. Thus, while 
T e g seems to imply the exploration of reversible jammed 
states within an energy ensemble, with 



P(E) 



Mr 



(18) 



describing the nature of the exploration, the validity of 
the energy ensemble to describe granular matter in the 
absence of energy conservation remains an open question. 
Here, we set the analogous Boltzmann constant for grains 
equal to unity for simplicity. 

Noting the drastic difference between the bath temper- 
ature and the effective temperature in a granular system, 
we are inspired towards a more careful analysis of the 
energy ensemble in slowly driven granular systems. The 
work of Edwards has promoted the concept of a volume 
ensemble, where the free volume per grain in a static 
granular system replaces the energy as the conserved 
quantity of the non-equilibrium system, at a particular 
volume fraction @, [35| • The basis for using the volume 
ensemble stems from the ability to conserve volume in a 
given packing and additivity of volume per grain. Fur- 
ther, it is possible to explore the configuration of states 
at a fixed volume, via experiment or simulation. The sta- 
tistical mechanics is then derived using methods similar 
to Boltzmann statistics for equilibrium systems. From 
these methods, one can obtain the compactivity, X, as 
a derivative of the entropy with respect to the volume, 
enabling the calculation of an equation of state in the 
volume ensemble as follows. 



X~ l = 



dS_ 
dV 



(19) 



The compactivity, X, is thereby assumed to be an 
equilibrium measure of a system within the framework 
of the volume ensemble, much like the bath tempera- 
ture of the energy ensemble. This assumption can be 
realized by performing an ABC experiment and testing 
a zero-th law of thermodynamics for volumes [4.0] . Ac- 
cording to the zero-th law of thermodynamics, if system 
A and C are in thermal equilibrium with system B re- 
spectively, then A and C are in thermal equilibrium with 



equal temperature. In granular system, such an exper- 
iment would require two granular systems with distinct 
volumes, Vi and V2, with the same X. Bringing these 
two systems together should result in a granular system 
of volume V = V\ + V2, at the same X, if the assump- 
tion is valid. This experiment is feasible due to the fact 
that it is always possible to prepare a system at a given 
volume fraction and will be the subject of future study 
and experiment, facilitated by recent theoretical findings 

Similar to the conservation of volume, boundary stress 
may also be a conserved quantity in jammed granular 
systems, and Edwards statistical mechanics for volume 
distributions could be applied analogously to the distri- 
bution of boundary stresses, II, or forces, referred to as 
the force ensemble [42|, |43[. The angoricity, A, is calcu- 
lated as the derivative of the entropy with respect to the 
boundary stress, and an additional equation of state is 
thereby achieved as follows: 



A ~ dU 



(20) 



This result can be combined with that of the volume 
ensemble in an effort to accurately define the statistical 
mechanics of static jammed granular matter. Such an 
approach remains a topic of ongoing research. 

However, slowly driven granular systems introduce yet 
another ensemble, the energy ensemble, from which the 
above defined T e g is derived. While the above results re- 
veal that T c g does not tend to zero as the magnitude of 
the driving force decreases, indicating extrapolation to a 
non-zero static quantity, it remains unclear how the ef- 
fective temperature may relate to the compactivity and 
angoricity as defined by Edwards statistics. Are we defin- 
ing a new static quantity by determining the static limit 
of T e ff, or are we expanding the statistical mechanics of 
jammed granular matter to include dynamic systems by 
relating T e g, X and Al T c g is obtained in the quasi-static 
limit 7 e — > + , while the volume and force ensembles cor- 
respond to 7 e = 0, exactly. Is it possible that a relation 
between T e g, X and A can be expected? 

There exists the further requirement of energy con- 
servation for the validity of a Boltzmann approach that 
would guarantee: 



TZ 



dS_ 
dE 



(21) 



As discussed above, energy is constantly dissipated in 
a driven granular system, through Coulomb friction and 
path-dependent tangential forces between grains. How- 
ever, the input of energy by the external driving force 
brings the system to a steady state where the average 
energy is constant over the time-scale of the experiment. 
This steady state energy could be likened to the con- 
served variable in a statistical formalism depicted in Eq. 
I21[ thereby introducing a thermodynamic meaning for 
T e s. 
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In a compressed emulsion system the absence of 
Coulomb friction and inter-particle tangential forces 
greatly simplifies the formalism *44|. Jamming occurs 
due to osmotic pressure, and the system remains athcr- 
mal as a result of the large particle size. A well defined 
potential energy exists due to the absence of tangential 
forces, corresponding to the deformation of the particles 
at the inter-particle contact points. Therefore, a restric- 
tion to use the energy ensemble in an effort to describe 
a jammed system is lifted, as frictional tangential forces 
no longer hinder energy conservation. 

Computer simulations of frictionless emulsion droplets 
[40( incorporate a simulated annealing method employ- 
ing an auxiliary temperature to sample the available 
jammed configurations. The simulated annealing method 
assumes a Boltzmann distribution, or a flat average as- 
sumption, over the jammed states of the emulsion. The 
Toff obtained by Eq. 1211 with simulated annealing meth- 
ods [J] is very close in value to the T a g obtained via the 
FDT calculations as in the present work. Such a result 
could indicate that ergodicity holds in this frictionless 
system, a further justification of the methods presented 
herein. Therefore, a firmer basis for the validity of us- 
ing the T s obtained in the quasi-static limit to describe 
the statistical mechanics of the same system at the static 
limit is achieved. 

Further, T e ff remains approximately constant with 
varying tracer particle size, implying a zero-th law of 
thermodynamics for slowly sheared jammed granular sys- 
tems. These statements further provoke the necessity 
for an ABC experiment, to test the zero-th law for 
the effective temperature, as well as similar experiments 
for the compacitivity and angoricity. Such experiments 
may enlighten us to understand under what conditions 

P(E) ~ e _, cff, P(V) ~ e"* and P(U) ~ e"§ may 
be valid in describing the statistics of the jammed and 
nearly jammed granular systems. 

At this point, we believe that the most prominent di- 
rection is the exploration of the volume and pressure en- 
sembles. Our understanding is that these ensembles may 
be sufficient to characterize the jammed state of granular 
matter, while the energy ensemble may be necessary for 
slowly moving granular systems. These are open ques- 
tions at the present time. Recent papers in the theory 
and simulation front suggest that the compactivity char- 
acterizes the system into a phase diagram at the isostatic 
point, while the angoricity will be necessary to describe 
the pressure ensemble of compressible granular matter 

El. 



V. SUMMARY 

In summary, this study focuses on the dynamics of 
slowly sheared granular matter in a 3D Couette cell. A 
mixture of spherical, transparent and bi-disperse grains 
are confined between two cylinders, having walls rough- 
ened by glued identical grains, with the inner cylinder 



rotated via motor. We compact the grains by means of 
an external pressure in the negative z-direction. Fluid 
matching the density and refractive index of the grains 
partially fills the cell, allowing tracking of tracer particle 
trajectories as a function of time. Tracers of varying den- 
sity and size are used. Multiple cameras track the tracer 
particle positions relative to the cylinders. 

We find that the angular velocity of the tracer parti- 
cles, u>e(r), follows an exponential relation with r, defined 
by the type of packing and geometry of the Couette cell. 
The velocity of the last layer of grains is non-zero, such 
that the shear band is located at the outer cylinder and 
ensures no formation of shear bands in the bulk. Near the 
outer cylinder LUg(r) decays slowly with increasing r, such 
that ug(r) can be approximated linearly with a constant 
local shear rate. The constant local shear rate ensures 
that the mobility and diffusivity of tracers, dependent 
on local shear rate, remain approximately independent 
of r. We define this region the "constant mobility and 
diffusivity region" , or the CMD region. 

An "effective temperature" , T c g , is realized by a 
fluctuation-dissipation relation generalized to granular 
materials. Statistical measurements are confined exclu- 
sively to the CMD region. The mobility in the vertical 
direction, M z , is found to be proportional to the shear 
rate, j e , for small enough values of %. As D z is also 
found proportional to shear rate, collapsing all the D z (r) 
for various shear rates shows a plateau in the CMD re- 
gion. An approximately constant effective temperature 
is obtained from measurements of the mobility and diffu- 
sivity, under a constant external applied force, and with 
sufficiently small shear rates. This effective temperature 
is calculated by an analogous equation used in equilib- 
rium statistical mechanics. We find this effective tem- 
perature to be independent of the tracer particle prop- 
erties, and dependent only on the packing density of the 
system. While this result describes an intensive property 
of the system, it remains an important future study to 
test the effective temperature against the laws of ther- 
modynamics. More specifically, a test of the zeroth-law 
of thermodynamics with respect to these non-equilibrium 
jammed systems could expand the scope of T c g beyond 
that of an intensive quantity of a particular system. A 
well defined effective temperature in the radial direction 
may exist, though its existence would require a constant 
external force applied in the radial direction. 

The probability distribution of the displacements in 
the radial direction, P(Ar), reveals exponential fluctua- 
tions. The analysis of the fluctuations reveals a power 
law, sub-diffusive, process, (Ar 2 ) ~ At a , with a less 
than unity. A similar analysis for fluctuations in the an- 
gular direction reveal a super-diffusive process, (A0 2 ) ~ 
At , with j3 greater than unity. Lastly, the probability 
distribution of the displacements in the vertical direc- 
tion are found have a Gaussian distribution such that 
(Az 2 ) ~ At. It is this linearity that defines vertical dis- 
placement as a diffusive process, and allows for the use 
of the Fluctuation-Dissipation relation to calculate the 
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diffusivity in the vertical direction. We further discover 
a linear relationship between angular displacement and 
the time between measurements A8 e ~ At, such that all 
mean square fluctuations can be defined in terms of A8 e 
for the small shear rates of our experiments. 

In the CMD region, the linear approximation of ug(r) 
proportional to approximately constant external shear 
rate, j e , allows for the collapsing of all tracer particle 
velocity curves via dividing ujg(r) by the shear rate. This 
collapse reveals a periodic shape with a small amplitude 
and periodic length roughly equal to the grain size. The 
effect is shown to be weaker in packings with smaller size 
grains. We further apply this remarkable scaling feature 
v z (r) and v r (r), achieving similar results. 

It is important to note that the effective temperature, 
defined in this study for small shear rates, does not re- 
main constant as the shear rate increases. While previous 
studies have discovered an an increasing effective temper- 
ature via simulations, we have measured diffusivity and 
mobility separately in an effort to calculate T c g through a 
fluctuation-dissipation relation. We find the diffusivity in 
the z-direction remains approximately constant through- 
out the range of shear rates used in this experiment, while 
the mobility in the z-direction approaches a plateau, ex- 
clusively increasing T e ff. Such an effect in the radial di- 
rection would be of great interest for future studies in 
sheared granular dynamics. 

The nearly constant value of T e g with respect to vary- 
ing tracer particle size indicates that a zero-th law of 
thermodynamics for slowly sheared jammed granular sys- 
tems could be valid and prompts one to perform an ABC 



experiment, fully testing the zero-th law for the effective 
temperature. 

As we work towards a more complete description of 
the statistical mechanics of jammed granular matter, we 
strive to incorporate the varied statistical ensembles into 
one fundamental picture. These ensembles include the 
energy ensemble, as described herein, along with the vol- 
ume and force ensembles, as proposed by Edwards. Such 
an incorporation may link static quantities of compaci- 
tivity and angoricity, describing volume and force ensem- 
bles, respectively, to the dynamic effective temperature 
presented in this study, derived from the energy ensem- 
ble. The exact nature of the relation between such quan- 
tities remains an open topic. Ultimately, these quantities 
will help to develop a thorough statistical description for 
jammed granular matter and reveal an equation of state. 
A deeper topic of concern is the formation of a clear def- 
inition of energy in jammed granular matter. Energy is 
not conserved in frictional systems and it remains open 
to debate as to how one would incorporate energy into 
the statistical mechanics. 

One possible approach to describe the energy of 
jammed systems is to consider the similarities between 
the inherent structure formalism of glasses and the ex- 
ploration of jammed states in granular matter, at least 
for the case of frictionless granular systems. Inherent 
structures probe a network of potential energy basins 
within an energy landscape. Such an approach toward 
the jammed states of granular matter may assist in un- 
derstanding exactly what is meant by energy within the 
framework of a non-equilibrium system. 
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